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We discuss techniques for probing the effects of a constant force acting on cold atoms using 
two configurations of a grating echo- type atom interferometer. Laser-cooled samples of 85 Rb with 
temperatures as low as 2.4 fjK have been achieved in a new experimental apparatus with a well- 
controlled magnetic environment. We demonstrate interferometer signal lifetimes approaching the 
transit time limit in this system (~ 270 ms), which is comparable to the timescale achieved by 
Raman interferometers. Using these long timescales, we experimentally investigate the influence of 
a homogeneous magnetic field gradient using two- and three-pulse interferometers, which enable us 
to sense changes in externally applied magnetic field gradients as small as ~ 4 x 10~ 5 G/cm. We 
also provide an improved theoretical description of signals generated by both interferometer config- 
urations that accurately models experimental results. With this theory, absolute measurements of 
B-gradients at the level of 3 x 10~ 4 G/cm are achieved. Finally, we contrast the suitability of the 
two- and three-pulse interferometers for precision measurements of the gravitational acceleration, g. 



I. INTRODUCTION 

Atom interferometers (AIs) have been employed to 
investigate a host of inertial effects over the past few 
decades. Such effects include the acceleration due to 
gravity [H-IH], gravity gradients 6-8], and rotations 8- 
ll| . Raman interferometric measurements of gravity 1 - 
3| use cold atoms and transit time limited experiments in 
an atomic fountain to reach a precision of ~ 3 parts per 
10 9 (ppb) with 1 minute of interrogation time. This tech- 
nique requires two phase-locked lasers to drive Raman 
transitions between two hyperfine ground states. It also 
requires state selection into the mp = magnetic sub- 
level to avoid sensitivity to £?-fields and _B-gradients, as 
well as velocity selection to guarantee that all interfering 
atoms have the same initial sub-recoil velocity. 

In contrast to the Raman interferometer, the grating 
echo- type AI [l2l . [l3[ uses a single off- resonant excitation 
frequency that drives a cycling transition with the same 
initial and final state. This AI requires no state or ve- 
locity selection, and is insensitive to both the AC Stark 
effect and the Zeeman effect [lij]. Additionally, as we 
will show on the basis of a theoretical model, the inten- 
sity of the AI signal is insensitive to uniform B-gradients 
provided the atoms are pumped into a single magnetic 
sub- level — which need not be mf = 0. 

In this work, we use two configurations of the grat- 
ing echo AI to demonstrate experiments with timescales 
comparable to those of Raman AIs. An improved the- 
oretical description of the echo AI has enabled accurate 
modeling of experimental data from which sensitive mea- 
surements of an externally applied B-gradient can be ex- 
tracted. This model is sufficiently general to describe 
all time-domain configurations of grating echo AIs, while 
accounting for a constant force on the atoms, as well as 
the sub-level structure of the atomic ground state. Re- 
cent work [15| with a particular configuration of this AI, 
has shown that measurements of the phase of the electric 
field are less sensitive to mirror vibrations. Investigations 



of the influence of £?-gradients on this AI configuration 
validate our predictions of gravitational effects, and in- 
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FIG. 1. (Color online) Recoil diagrams for the two-pulse (a) 
and three-pulse (b) AIs in the absence of any external forces. 
Standing wave excitations are labeled by SWl - SW3. At 
t = Ti the atom is diffracted into a superposition of momen- 
tum states differing by integer (n) multiples of hq. A second 
sw pulse is applied at t — T2 which further diffracts the atoms. 
In the two-pulse case, interference between states differing by 
hq (An = ±1) occurs at times t^J ho = Ti + (N+l)T 2 i. For the 
three-pulse scheme, a third sw pulse is applied at t = T3 which 
produces interference at times = Ti + (N + l)T2i + Ts2- 
Examples of interfering trajectories for N = 1 and 2 are la- 
beled by solid black and dashed blue lines, respectively. Cir- 
cles indicate locations where interference fringes occur with 
spatial frequency q. 
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dicate that this AI is particularly well-suited for precise 
measurements of the gravitational acceleration, g. 

We begin with a review of the two AI configurations 
used in this work, which are illustrated in Fig. [T] The 
two-pulse echo AI [12, 0, [l|| [ljj utilizes short (Raman- 
Nath) standing wave (sw) pulses to diffract a sample 
of laser cooled atoms at t = T\ into a superposition 
of momentum states: \nhq). Here, n is an integer and 
q = ki — k 2 ~ 2fc is the difference between the traveling 
wave vectors comprising the sw. At t — T 2l a second sw 
pulse further diffracts the atomic wave packets — creating 
sets of center-of-mass trajectories that overlap and pro- 
duce interference in the form of a density modulation in 
the vicinity of tgL = T! + (iV+l)T 21 , where T 21 = T 2 -T x 
and N — 1,2,... is the order of the echo, as shown in 
Fig. 1 (a) The induced density modulation is coherent for 
r co h = 2/qa v ~ 3 fis about these "echo" times, beyond 
which the modulation dephases due to the distribution of 
velocities in the sample. Here, a v = (2/csT/M) 1//2 char- 
acterizes the width of the velocity distribution along z. 

A traveling wave pulse is applied along the ^-direction 

(2) 

in the vicinity of L to "read out" the amplitude of 
the grating by coherently Bragg scattering light along 
the — ^-direction. The duration of this signal is limited 
by the coherence time, r co h. Due to the nature of Bragg 
diffraction, this back-scattered light is proportional to the 
Fourier component of the density distribution with spa- 
tial frequency q. This harmonic is only produced by inter- 
ference of momentum states that differ by hq (An = ±1). 
As a result, the two-pulse AI exhibits a temporal modu- 
lation at the atomic recoil frequency, uj q — hq 2 /2M, and 
is therefore sensitive to recoil effects. 

The three-pulse "stimulated" grating echo AI (hence- 
forth referred to as the three-pulse AI) was first demon- 
strated in Ref . , 13 using a single hyperfine ground state, 
and was termed a "stimulated" echo due to similari- 
ties in pulse geometry with the stimulated photon echo 
scheme 18-21]. Recent work involving this interferom- 
eter [ll| has shown certain advantages over the two- 
pulse scheme for phase measurements of the atomic grat- 
ing. The three-pulse AI involves applying two sw pulses 
at t = Ti and t = T2, followed by a third pulse ap- 
plied at t = T 3 = T 2 + T32, where T 32 = T 3 - T 2 . 
This pulse geometry produces an echo in the vicinity of 
€io — Ti + (N + 1)T 21 + T 32 , as shown in Fig. [lfb) 



However, unlike the two-pulse AI where all pairs of tra- 
jectories produced by the second pulse interfere at the 
echo times, for the three-pulse AI only momentum states 
of the same order (An = 0) after the second pulse pro- 
duce interference at the echo times for arbitrary T 2 i and 
T32. For this reason, the signal produced by this interfer- 
ometer as a function of T32 (with T 2 \ fixed) is insensitive 
to atomic recoil effects (i.e. no temporal modulation) and 
is therefore ideal for probing other effects — such as those 
due to a constant force on the atoms. 

Reference H3 extensively reviews the grating echo AI 
and discusses applications relating to atomic recoil [lg, 
UtI |23T - [26| , gravity and magnetic gradients (2?j] ■ 



Previous experiments based on this AI [l2j, [l3|, [lfj, Il7l . 
[H-[H [H were typically limited to T 21 < 10 ms by de- 
coherence effects due to spatially and temporally varying 
B-fields. Additionally, the sample temperature (typically 
~ 50 uK) and excitation beam configuration (fixed fre- 
quency sw with ~ 0.5 cm diameter) limited the transit 
time in these experiments. In this work, we have im- 
proved the level of B-field and _B-gradient suppression 
by using a non-magnetic vacuum chamber, which has 
enabled the extension of AI signal lifetimes. The mag- 
netically controlled environment allows a sample of 85 Rb 
atoms to be cooled to temperatures as low as 2.4 uK. By 
expanding the excitation beam diameter to ~ 2 cm, and 
chirping the sw pulses to cancel Doppler shifts, echo AI 
signal lifetimes of ~ 220 ms and transit times of ~ 270 
ms have been achieved. These timescales are comparable 
to those of fountain experiments involving Raman AIs 
[E EL [HI- In contrast, long-lived echo AI signals have 
been observed only by using magnetic guides to limit 
transverse cloud expansion [15j. 

The experimental apparatus presented here has made 
it possible to exploit the aforementioned advantages of 
the echo AI for a variety of precision measurements, such 
as the atomic recoil frequency [26| and the gravitational 
acceleration [22| . that are currently underway. Addi- 
tionally, we recently utilized this apparatus to perform 
a coherent transient experiment with cold Rb atoms to 
achieve a precise determination of the atomic g-factor 
ratio [30l |. 

In this Article, we apply long timescales to understand- 
ing and detecting the effects of S-gradients using the 
two-pulse grating echo AI, as well as a three-pulse "stim- 
ulated" grating echo AI [l3|, [l5| . The passive detection 
of magnetic anomalies is of interest for various applica- 
tions, such as submarine and mine detection where the 
ambient magnetic noise of the environment is large com- 
pared to the sensitivity of the instrument [3l| . The influ- 
ence of gravity and B-gradients on AI experiments has 
been considered in the past. Reference [H calculates how 
such forces affect the visibility of interference patterns in 
atomic diffraction experiments. In previous work (2?| . we 
demonstrated the effect of both gravity and i?-gradients 
on the two-pulse AI. A theoretical description of these 
effects based on a spin- 1/2 system was able to explain 
the basic signal dependence on the pulse separation, T 2 \, 
but was insufficient to model experimental data. 

This work relies on an improved theoretical description 
of a generalized echo AI that includes an arbitrary num- 
ber of sw excitation pulses, the effects of a constant force 
on the atoms, spontaneous emission and the sub-level 
structure of the atomic ground state (the 5Si/ 2 F = 3 
state of 85 Rb is used in the experiment). Coupled with 
these theoretical predictions, we achieve sensitivity to 
changes in B-gradients at the level of ~ 0.04 mG/cm. In 
addition, absolute measurements of £?-gradients as small 
as ~ 0.3 mG/cm, and sensitivity to the curvature of in- 
fields are demonstrated. These results are consistent with 
independent measurements of the spatial variation in the 
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£?-field using a flux-gate magnetometer. These studies 
help place limits on the sensitivity of a broad class of 
time-domain AIs to B-gradients. 

We also consider implications for achieving precise 
measurements of g using the two- and three-pulse echo 
AIs. In particular, analysis of the three-pulse AI suggests 
there are significant advantages for measuring g over the 
two-pulse AI. Although the experimental apparatus used 
in this work is not designed to detect gravitational effects, 
predictions of the grating phase modulation due to grav- 
ity for both AIs have been validated by measuring the 
effects of externally applied B-gradients. Measurements 
of g using these AIs will be presented elsewhere. 

This Article is organized as follows. In Sec. UU we 
present theoretical predictions for the two- and three- 
pulse AI signals in the presence of a homogeneous B- 
gradient. In Sec. IIII1 we describe details related to the 
experimental apparatus. We present the results of exper- 
iments related to long timescales in Sec. IIV| and discuss 
measurements of i3-gradients using both the two-pulse 
and three-pulse techniques. Section [V] discusses the fea- 
sibility of a precise measurement of g using the formal- 
ism developed to describe B-gradients. We conclude in 
Sec. EH The Appendix presents a calculation of the sig- 
nal generated by a generalized echo AI — encompassing 
the two- and three-pulse AIs — in the presence of a con- 
stant force. 



II. THEORY 

In this section, we present the key results of calcula- 
tions for both the two- and three-pulse AI signals in the 
presence of a homogeneous B- gradient. Details of the 
calculations — which are sufficiently general to account for 
any constant force on the atoms, and an arbitrary num- 
ber of excitation pulses — are presented in the Appendix. 

In general, the sensitivity of these interferometers can 
be characterized by the space-time area they enclose. 
Since only those states differing by hq at the echo time 
contribute to the signal, the area of both AIs is primarily 
controlled by T 2 \. In the absence of any external forces, 
the areas of the two- and three-pulse AIs can becalcu- 
latcd by inspecting their recoil diagrams [Figs. 1(a) and 
1(b) respectively] 
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where M is the mass of the atom. Henceforth, quan- 
tities containing superscripts (2) or (3) indicate the in- 
terferometer for which that quantity applies. At first 
glance, it might appear that the three-pulse AI encloses 
a larger area than the two-pulse AI due to the extra 
term in Eq. (Ilb|) . However, one must compare the en- 
closed areas at the same echo times, which are given by 
( lo = T 1 + (N+l)T^andt^ 



for the two- and three-pulse schemes, respectively. By 

setting t^L =jicL> it; can be shown that - = 
HqNT^ 2 /2M(N+l). This suggests that the two-pulse AI 
is always more sensitive to external forces than the three- 
pulse AI. Nevertheless, the three-pulse AI offers a unique 
feature: the spatial separation between interfering wave 
packets remains constant between the application of the 
second and third sw pulses. This is advantageous because 
larger spatial separations leads to increased decoherence, 
and therefore reduced timescale in the experiment [l5j . 
Since the separation can be precisely controlled by the 
pulse separation T21, one can increase the signal lifetime 
by using smaller T%\. 

Additionally, since the signal generated by the two- 
pulse AI is modulated at the recoil frequency, u g , there 
are periodic regions where the signal-to-noise ratio is less 
than one and not well-suited for accurate phase measure- 
ments. However, the three-pulse technique is insensitive 
to atomic recoil if T21 is fixed. Therefore, the scattered 
field amplitude has no additional modulation at co q as T32 
is varied — allowing regions of low signal-to-noise ratio to 
be avoided. 

Both the gravitational force and a constant B-gradicnt 
produce a constant force on the atoms, J- = !Fz, which 
generates a phase shift in the atomic interference pattern. 
The basic physical mechanism that produces this phase 
shift is a difference in potential energy between the two 
arms of the AI. One can compute the relative phase be- 
tween the two arms Acf> = (6>b — <Sa)/^ using the classical 
action Q 



s(t)= few), 

Jo 



(2) 



where C = Mz 2 /2 + Tz is the Lagrangian in this case. 
If 6>b and 6>a represent, respectively, the action along 
the upper and lower arms of the two-pulse AI, it can be 
shown that the phase shift between these arms is 



A0 (2) = N(N + 1) 



M 



-T, 



21 



+ N 2 qv T 21 , (3) 



where Vq is the initial velocity of the atom along the 
^-direction. The term proportional to vq is due to the 
relative Doppler shift between the two arms of the AI. 
Since the atomic sample has a finite velocity distribution 
(characterized by a 1/e radius, a v , and temperature, T), 
this term is responsible for the coherence time of the echo: 
T co h = 2/qa v . As expected, the contribution to the phase 
shift from the potential energy (the term proportional to 
J 7 ) is independent of the initial velocity of the cloud. 

A similar calculation for the relative phase shift be- 
tween the arms of the three-pulse AI yields 



t 



Ti + (AT+1)T 2 ^ andi^ = r 1 + (A+l)T 2 ( ? ) +T 32 
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+ Nj^T 32 T 21 + Nqv (T 32 + NT 21 ) 
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This expression is similar to Eq. Q, with additional 
terms proportional to the pulse spacing T 32 . One can 
vary either T 2 i or T32 to detect phase modulation pro- 
duced by an external force, F . However, since there are 
no terms containing the phase uj q T 32l one can effectively 
turn off the sensitivity to atomic recoil by fixing T21. 
This makes the three-pulse AI ideal for investigating the 
effects due to F, especially when qFT 21 /M >• uo q since 
no additional modulation at w q is present. This is par- 
ticularly advantageous for measurements of gravity, as 
discussed in Sec. [V] 

Since the two-pulse AI is intrinsically more sensitive 
than the three-pulse AI, it is better suited to measure- 
ments of F when qFT 2 \/M < Lu q . In this work, we 
demonstrate this feature by measuring externally applied 
B-gradients. The sensitivity of the three-pulse AI to F 
can be enhanced by utilizing the additional phase pro- 
portional T32T21 in Eq. Experimentally, this can be 
accomplished by varying both pulse separations, T21 and 
T32, with T21 varied in integer multiples of the recoil pe- 
riod: T q = 1T/U! q . 

To determine the response of the grating echo AI in 
the presence of a constant force, we assume a potential 
energy with the form 



U{z) 



-Mz, 



(5) 



where Ai = —dU/dz is a matrix operator with units 
of force that commutes with both the position (z) and 
momentum (p) operators. In the case of a constant B- 
gradient, the potential is 



U(z) 



-H-B(z) 



F 7 z, 



(6) 



where gp is the Lande g-factor, [Ib is the Bohr magne- 
ton, B{z) — (3z is the magnetic field vector with gradient 
P along the z-direction (also assumed to be the quanti- 
zation axis), and F z is the projection operator for total 
angular momentum, F. In this case, A4 — FF z /h and 
the force is F = gp^sP, where F z operates on the basis 
states \Fmp) and has eigenvalues hmp. 

In both interferometer schemes, the phase of the grat- 
ing is imprinted on the electric field back-scattered from 
a traveling- wave read-out pulse applied in the vicinity of 
an echo time. For the two-pulse AI, measuring the phase 
of the scattered field is equivalent to measuring the rela- 
tive position of the grating, since both the position and 
the phase scale as T 21 . Similarly, in the three-pulse case, 
the phase measured as a function of T32 is proportional 
to the velocity of the grating — which scales as T21. 

We first examine the effects of -B-gradients on the two- 
pulse AI, followed by a comparison with the three-pulse 
AI. 



A. Two-Pulse Interferometer 

In general, the electric field scattered by the atoms 
at the time of an echo is proportional to the amplitude 



of the Fourier harmonic of the atomic density grating 
with spatial frequency q. For the case of an external B- 
gradient, (3, the scattered field has distinct contributions 
from each magnetic sub-level: 



Ef (t; T)=J2 E%1 (*; T)e im ^ < t;T 



(7) 



( 2) 

where Em F is the field scattered by the state \Fmp) 
[given by Eq. (|A21[) ] and mp(f>^ is the phase shift of 
the density grating produced by the same state in the 
presence of the B-gradient. For the ./V th order echo at 
*ecL = T i + (N + 1 ) I 2i ) with the set of onset times 
T = {Tx,Tx + T 21 } and At = t - t^ , the phase shift of 

(2) 

the grating tpi '(t; T) is given by 



^ 2) (Ai;T 21 ) 



{N(N + i)Tl 



2[F + (N 



2M 

Hl)T 21 ]Ai + Ai 2 } 



(8) 



The general form of this equation for a constant force, F, 

is given by Eq. (|A23|) in the Appendix. In the discussions 

that follow, we take At = which corresponds to the echo 

time. Since the echo lasts for r co h ~ 3 /Lts about At = 0, 

the signal is obtained by integrating the back-scattered 

field over this time. 

Equations ([7]) and (J5|) indicate that the field amplitude 

scattered from state \F mp) exhibits phase modulation 

(2) 

as a function of T21 at a frequency mpujp (T21) due to 
the presence of the gradient, where 



dT 21 



M 



(9) 



JV(JV + l)T 3 i + (N+ l)At] 



This modulation frequency depends linearly on /3 and the 
pulse separation, T 2 i (i.e. the frequency is chirped with 
T21). The phase modulation of the grating produced by 
state \F mp) also scales linearly with the magnetic quan- 
tum number, mp, as shown in Eq. ([?])• For an arbitrary 
set of magnetic sub-level populations, the total scattered 
field [Eq. j7|] contains all harmonics mpuip {T 2 i), where 
mp = —F, . . . , F. If more than one sub-level is popu- 
lated, interference between the fields scattered off of each 
state produces modulation in the total scattered field. 
This effect can then be detected in the field amplitude, 
Ep, or the field intensity, l-E^ 2 '! 2 , by varying (3 or the 
pulse separation, T21. The amplitude of each harmonic 
comprising this modulation is determined by the sub- 
level populations, as well as the transition probabilities 
between ground and excited state sub-levels. 

If the system is optically pumped into a single sub- 
level, such as the extreme state: \F F), then the phase 
modulation of the grating only affects the phase of the 
electric field — which cannot be observed using intensity 
detection. Instead, one can use heterodyne detection to 
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FIG. 2. (Color online) The predicted B-gradient signal for 
the two-pulse AI (a,b) [based on Eq. Q], where T21 is varied 
in steps of r q ~ 32 /is, and the three-pulse AI (c,d), where 
T32 is varied. In all plots the solid red line is the real part of 
the scattered field, the black dashed line shows the field inten- 
sity (which is measured in the experiment), the B-gradient is 
fixed at P = 10 mG/cm, and the first order echo (TV = 1) is 
used. In parts (a,c) the excitation beams are circularly po- 
larized (|(/l| = 1) and the magnetic sub-level populations are 
equally distributed among the F — 3 ground state of 85 Rb. 
In this case, the field undergoes amplitude modulation with 
multiple frequency components, and the field intensity ex- 
hibits modulation with a contrast < 100%. For the two- pulse 
AI (a), the field intensity exhibits modulation at a chirped 
frequency, whereas for the three-pulse AI (c) the modulation 
is not chirped. Parts (b,d) show the field generated by a 
sample that is optically pumped equally into the |3 —3) and 
|3 3) states with linearly polarized excitation beams (§l = 0). 
Here, both the field and the field intensity exhibit modulation 
with only one frequency component, and the contrast of the 
oscillations is 100%. For both (c,d), T21 was fixed at a typical 
experimental value of 5 ms. 



measure the electric field amplitude and obtain the rel- 
ative phase of the scattered light [HI, [2^, [27| ■ Further- 
more, if the system is optically pumped into the \F 0) 
state, there is no phase modulation due to £?-gradients 
since this state is insensitive to magnetic fields. 

Figures 2(a) and |2(b)] show the expected two- pulse AI 
signal as a function of T21 in steps of the recoil period, 
T g = ix/uiq (~ 32 /is for 85 Rb). Since lo^ < uj qi incre- 
menting T21 in this fashion eliminates additional modula- 
tion due to atomic recoil. Figure 2(a)| shows the signal for 
equally distributed sub- level populations, while Fig. 2(b) 



is for an optically pumped system in the two extreme 
states: |3 —3) and |3 3). Both of these figures show am- 
plitude modulation, but in the optically pumped case 
there is only one frequency component present and the 
modulation occurs with maximum contrast — increasing 
the sensitivity to gradients. 

Eliminating the amplitude modulation in the sig- 
nal due to S-gradients [shown by the dashed lines in 
Fig. 2(a) is a key requirement for precision measure- 



ditions can be realized with sufficient suppression of am- 
bient i3-gradients in a glass cell. It is also possible to 
eliminate sensitivity to i3-gradients using intensity de- 
tection if the atoms are pumped into a single magnetic 
sub-level. 



B. Three-Pulse Interferometer 

The effects due to B-gradients manifest themselves dif- 
ferently in the three-pulse interferometer. We derive the 
expression for the signal in the Appendix [see Eqs. (|A25|) 
and (j A29|) ] and find that the amplitude of the scattered 
field does not depend on the time between the second 
and third sw pulses, T32, but only on T21 — similar to the 
two-pulse AI. However, the phase of the grating in the 
three-pulse case depends on both T 2 \ and T32: 

q9F ^{N(N 



bf\At; 



+ 2[Ti 

In this case, 
T = {Ti,Ti4 



T 
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l)Tai]At + At 2 } 



(10) 



the set of pulse onset times is given by 
T 21 ,T 1 +T 21 +T 32 } and At = t-t^ ho . This 
phase is identical to Eq. for the two-pulse interferom- 
eter with the addition of the two terms proportional to 
T 32 . Equation (fit))) suggests that the force can be deter- 
mined by measuring the phase modulation of the grating 
as a function of either T21 or T32, or by varying both 
pulse separations simultaneously. Varying T32 produces 
a phase modulation of the atomic grating at a frequency 
that is proportional to T 2 \\ 
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(T21) 
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ments of u>„. We will show in Sec. IIVI that these con- 



Figures 2(c) and |2(d)| show the expected three-pulse 
signal as a function of T32, with T 2 \ fixed at a typical ex- 
perimental value of 5 ms, in the presence of a S-gradient 
P = 10 mG/cm. When the sub-level populations are 
equally distributed [Fig. |2(c)] the phase of the total scat- 
tered field contains multiple frequency components — one 

(3) 

for each sub- level: mpui^ . The interference between 
these components produces a modulation in the total 
scattered field amplitude. This is similar to the two-pulse 
case shown in Fig. 2(a)| except that the modulation oc- 
curs at a single frequency that is fixed by /?, T 2 \ and 
N . For a sample that is optically pumped equally into 
the two extreme states: |3 —3) and |3 3), as shown in 
Fig. |2(d)| there is only one frequency component present 
in the scattered field. In this case, the amplitude mod- 
ulation occurs with greater contrast than for any other 
configuration of sub-level populations. 



III. EXPERIMENTAL SETUP 

We now review the experimental setup that has made 
possible long-lived grating echo AI signals. This setup 
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is substantially different from previous echo experiments 
[IH HH HH l27| after implementing many improvements. 
These include suppression of stray magnetic gradients 
using a non-magnetic chamber, increasing the trapped 
atom number with large diameter beams, extending the 
transit time by cooling the sample to ~ 10 /iK and imple- 
menting large diameter excitation beams, and by chirp- 
ing the excitation frequencies to eliminate Doppler shifts 
associated with the falling cloud. 

The experiment utilizes a sample of laser-cooled 85 Rb 
atoms in a magneto-optical trap (MOT) containing ap- 
proximately 10 9 atoms in a Gaussian spatial distribution 
with a horizontal e" 1 radius of ~ 1.7 mm. The MOT 
is contained in a borosilicate glass cell maintained at a 
pressure of ~ 1 x 10~ 9 Torr. In addition to the anti- 
Hclmholtz coils used for trapping, three pairs of square 
quadrupole coils are centered on the MOT, as shown in 
Fig-El Each square frame contains two overlapping coils, 
one connected in the Helmholtz configuration with the 
coil in the opposite frame, and the other in the anti- 
Hclmholtz configuration. These sets of coils are used to 
cancel ambient magnetic fields and field gradients over 
the volume of the MOT at the level of ~ 1 mG and 
~ 0.1 mG/cm, respectively. The initial set points for the 
currents in the canceling coils that produced ~ 1 mG 
of B-field suppression were determined using an atomic 
magnetometer experiment (30j that allowed the field at 
the location of the MOT to be measured. 

Light derived from a Tksapph laser with frequency vl 
and linewidth ~ 1 MHz is locked to the 5S]/ 2 F = 3 — s- 
F' = 4 transition (vl — vq) using saturated absorption 
spectroscopy. The light is then shifted 130 MHz above 
resonance by an acousto-optic modulator (AOM) operat- 
ing in dual-pass mode such that vl = v$ + 130 MHz. A 
separate "trapping" AOM shifts this light by -148 MHz 
such that the detuning is A = —14 MHz (vl = — 14 
MHz). Approximately 370 mW of this light is transmit- 
ted through an anti-reflection-coated, single-mode optical 
fiber (operating at 60% efficiency) and expanded to a di- 
ameter of ~ 5.4 cm for trapping atoms from background 
vapor. 

An external cavity diode laser is used to derive repump 
light for the trapping setup. It is locked to the 5Si/ 2 
F = 2 — > F' = (2, 3) crossover transition and up-shifted 
by ~ 32 MHz using an AOM. Approximately 25 mW of 
repump light is obtained after coupling through the same 
optical fiber as the trapping light. At t = 0, the MOT 
coils are pulsed off in ~ 100 /is, while the trapping and 
repump beams are left on for 6 ms of molasses cooling. 
For ~ 3 ms of this time, the detuning of the trapping light 
is linearly chirped from A = — 14 MHz to —50 MHz to 
further cool the atoms, and the power is simultaneously 
ramped down in order to reduce heating due to photon 
scattering. With this procedure we achieve temperatures 
as low as T — 2.4 /xK. 

Light from the Ti:sapph laser is also used to derive the 
AI pulses. A "gate" AOM operating in dual-pass con- 
figuration shifts the undiffracted light from the "trap- 




FIG. 3. (Color online) Diagram of the experiment. The ex- 
citation beams are both cr + -polarized by the A/4-plates. The 
glass cell has approximate dimensions 7.6 x 7.6 x 84 cm. Each 
pair of square quadrupole coils have a side length of ~ 66 
cm and contain overlapped coils wired in both Helmholtz and 
anti-Helmholtz configurations for canceling B-fields and B- 
gradients, respectively. 



ping" AOM from v L = v Q + 130 MHz to v L = u + 290 
MHz. The gate AOM is also pulsed so as to serve as 
a high-speed shutter during the experiment. The light 
from the gate AOM is split and sent into two separate 
AOMs (referred to as the "fei" and "fc 2 " AOMs) operat- 
ing at 240 MHz ± 5(t) that produce the sw pulses. Here, 
S(t) = gt/X is a time-dependent frequency shift that is 
added to (subtracted from) the radio frequency (rf) driv- 
ing the k\ (fc 2 ) AOM using an arbitrary waveform gener- 
ator, as shown in Fig. [JJ Chirping the excitation pulses 
in this manner cancels the Doppler shift of the atoms 
falling under gravity. The rf driving these AOMs is also 
phase locked to a 10 MHz rubidium clock to eliminate any 
electronically induced phase shifts. Light entering the k\ 
AOM is downshifted by 240 MHz - 5(t) and sent into an 
optical fiber that carries the light toward the MOT. Simi- 
larly, the k 2 AOM downshifts the light by 240 MRz+5(t). 
In this configuration, the detuning of the fci (fc 2 ) pulse 
is Ai = 50 MHz - S(t) [A 2 = 50 MHz + £(*)]. This light 
is coupled into a separate fiber and aligned through the 
MOT along the vertical direction, as illustrated in Fig. fj] 
The output of both fibers is expanded to a e~ 2 diameter 
of ~ 2 cm. The rf pulses driving the k\ and fe 2 AOMs 
are controlled using TTL switches with an isolation ratio 
of 100 dB, which produces optical pulses with rise times 
of ~ 20 ns. The "gate" AOM is turned off between ex- 
citation pulses to further reduce background light from 
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FIG. 4. (Color online) Schematic of the rf chain used for 
chirped AI pulses. A phase-locked loop (PLL) generates a 
220 MHz rf signal which is split and mixed with the output 
of two separate arbitrary waveform generators (AWGs). The 
AWGs are triggered at the start of the experiment to output 
a frequency sweep from 20 MHz to 20 MHz±<5(t) after a time 
t, where S(t) — gt/\, g ~ 9.8 m/s 2 and A ~ 780 nm. The sum 
frequency from the mixers is isolated using a band-pass filter 
(BPF) with a center frequency of 240 MHz and a 5 dB pass- 
band of 4 MHz. The outputs of the BPFs are pulsed using 
a set of transistor-transistor logic (TTL) switches, which are 
then sent to the fci and h,2 AOMs. The two-pulse AI sequence 
for both fci and Ai2 are shown. Here, PI and P2 refer to 
traveling wave components comprising the first and second sw 
pulses, and RO denotes the traveling wave read-out pulse sent 
along fci . Both AWGs and the PLL are externally referenced 
to a 10 MHz Rb clock. 



reaching the atoms. 

In the vicinity of any given echo (see Fig.Q]), the read- 
out pulse is applied to the sample along the fci-direction 
and a coherent back-scattered field from the atoms occurs 
along the direction of . The power of the scattered field 
is recorded as a function of time using a photo-multiplier 
tube (PMT) that is gated on for 9 /is. The echo signal 
lasts r co h ~ 3 fis before coherence is lost due to Doppler 
dephasing. For < 10 ms, the scattered field can reach 
powers greater than 100 /iW. However, for T21 > 10 ms, 
the signal size decreases exponentially. The noise floor 
for the PMT is approximately 0.1 /xW. Typically, one 
computes the time-integrated area of the echo signal as 
a measure of the signal size for a given set of parameters. 
Since this quantity has units of energy, it is henceforth 
referred to as the echo energy. 



IV. RESULTS AND DISCUSSION 

We now review the main experimental results of this 
work relating to long-lived AI signals and sensing exter- 
nally applied £?-gradients. 



A. Investigations of AI Timescale 

Figure Eta) shows a measurement of the temperature 
of the laser cooled sample. At t = To, all optical and 
magnetic fields associated with the MOT are switched 
off and the atoms are allowed to thermally expand in the 



dark. At t = Tq + T cxp , the trapping and repump beams 
are turned back on and a calibrated charged-coupled de- 
vice (CCD) is triggered to photograph the cloud with an 
exposure time of 100 fis. This process is repeated for 
various expansion times, T oxp , and the e -1 radius of the 
cloud, i?, is measured by fitting to the Gaussian intensity 
profiles obtained from each image. The temperature is 
obtained by fitting to a hyperbola [HI, H3| with the form 
R = [i?o+erJ;(T CX p — to) 2 ] 1 / 2 , where R is the initial cloud 
radius, a v = (2fcsT/M) 1 / 2 is the e" 1 radius of the veloc- 
ity distribution and to is a phenomenological offset from 
7c X p = 0. The data shown in Fig. EJa) give a tempera- 
ture of T - 2.4 fiK in 85 Rb. This relatively low MOT 
temperature is attributed to the well-controlled magnetic 
environment within the glass cell, as well as the molasses 
cooling procedure described above. 

Measurements of the AI signal lifetime under different 
pulse configurations are shown in Fig. Elk), with each 
configuration explained schematically in Fig. Etc). For 
the transit time measurement, the two- pulse AI configu- 
ration was used with T%\ fixed. The excitation and read- 
out pulses were incremented synchronously. The signal 
lifetime for the three-pulse AI was determined by fixing 
T21 and varying the third sw pulse and read-out syn- 
chronously. For the two-pulse AI, the lifetime was mea- 
sured by fixing the first sw pulse and incrementing the 
second sw pulse and read-out in steps nT q and 2nr q , re- 
spectively, where n = 10. 

Here, all ambient S-fields and B-gradients are can- 
celed along all three axes at the level of ~ 1 mG and 
~ 0.1 mG/cm, respectively. The transit time data was 
obtained by using the two-pulse AI with T21 fixed at 
~ 1.690 ms and varying the time of all sw pulses relative 
to the time of trap turn-off, To- In this measurement, the 
AI signal is proportional to the number of atoms that re- 
main in the volume defined by the ~ 2 cm diameter exci- 
tation beams during the thermal expansion of the cloud. 
Although the echo energy spans almost three orders of 
magnitude as it decays exponentially, signals are clearly 
distinguishable from the noise floor (~ 0.1 pJ) at times 
as large as ~ 270 ms, as shown in Fig. Elk). This time 
represents the transit time limit for the conditions of our 
experiment — corresponding to a drop height of ~ 36 cm. 
This distance nearly coincides with the bottom viewport 
of the vacuum system. We emphasize here that such life- 
times are not possible with this interferometer unless the 
frequencies of the k\ and k,2 beams are oppositely chirped 
such that the Doppler shift due to gravity [S(t) — gt/X] 
is canceled or the bandwidth of the sw pulses is large 
enough to account for such a shift. The frequency chirp 
puts the sw pulses on resonance for the two-photon tran- 
sition back to the same ground state for all times during 
the sample's free-fall. 

The signal lifetime for the two-pulse AI configuration 
is shown as the red curve in Fig. Efb). Here, the signal 
lasts approximately 130 ms, corresponding to T21 ~ 65 
ms. To the best of our knowledge, this is the largest 
timescale observed with this interferometer, correspond- 
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FIG. 5. (Color online) (a) Temperature measurement of the laser-cooled sample. The horizontal e _1 cloud radius is measured 
as a function of expansion time using a CCD camera. A hyperbolic fit to these data (shown as the solid red line) yields a 
measurement of T = 2.4 ± 0.2 fjK. (b) Data showing the signal lifetime for various pulse configurations. The horizontal axis is 
the time of the read-out pulse, Tro, relative to the time of trap turn-off, To, which signifies the start of the experiment. An 
echo energy of ~ 0.1 pJ is equivalent to the level of noise in the detector. These data were obtained with pulse durations for 
the two-pulse (three-pulse) AI: Ti = 800 (800) ns, T2 = 100 (70) ns, T3 = (70) ns; pulse intensity / ~ 64 mW/cm 2 ; and a 
sample temperature of T ~ 20 [iK. Pulse separations for each configuration are shown in the figure, (c) Pulse configurations 
for the data shown in (b). The red arrows distinguish which pulses were varied and the time step, dT, indicates the amount 
each pulse was incremented relative to the others. For the two-pulse AI, dT = nr q , where n = 10. Pulses without arrows are 
fixed in time. 



ing to more than a factor of 6 improvement over our 
previous work [lj| [24], [2^, [27j. However, the lifetime of 
the two-pulse echo is still limited by decoherence from a 
small, inhomogeneous £?-gradient that the atoms sample 
over the ~ 8 cm they have fallen in 130 ms. A non-linear 
B(z) produces a spatially-dependent force between inter- 
fering trajectories — resulting in a differential phase shift 
between paths of the interferometer that causes dephas- 
ing and, therefore, a loss of signal. Such a non-linearity in 
B(z) has been measured to be d 2 B/dz 2 ~ —0.4 mG/cm 2 
with a flux-gate magnetometer placed at different spatial 
locations around the glass chamber. This curvature is 
produced by a combination of non-ideal coil configura- 
tions and the presence of nearby ferromagnetic materials. 

There are two important features that should be rec- 
ognized from the data for the three-pulse AI shown in 
Fig. E^b). First, at T RO - T ss 0, the echo energy for 
the three-pulse AI is a factor of ~ 2 smaller than that 
of the two-pulse AI. This comes about because the ad- 
ditional Kapitza-Dirac pulse involved in the three-pulse 
AI produces fewer pathways that result in interference 
at the echo time compared to the two-pulse AI. Second, 
the lifetime of the three-pulse echo depends strongly on 
the value of T21. As T21 increases, the signal lifetime 



approaches that of the two-pulse AI. This feature comes 
about because, between the second and third sw pulses, 
the wave packets that interfere at the echo times have a 
constant spatial separation [see Fig. 1(b)] , which is given 
by Az = NhqT2i/M . From this expression, it is clear 
that Az can be controlled by T21 and the choice of echo 
order, N. By decreasing this separation, the interferome- 
ter becomes less sensitive to decoherence from non-linear 
B-fields since phase shifts produced by this effect be- 
come approximately common mode between interfering 
momentum states. Reference [TBI also used this interfer- 
ometer and a magnetic guide to show that smaller spatial 
separations lead to increased timescales. 

In general, the lifetime for the three-pulse echo can 
be tailored to last much longer than that of the two- 
pulse echo, which is advantageous for precisely measuring 
the effects of external forces. For example, we achieve 
timescales as large as ~ 220 ms for T21 fixed at ~ 1.3 
ms — which is much closer to the transit time limit than 
the lifetime of the two-pulse echo. To the best of our 
knowledge, the only experiment that has achieved longer 
timescales for this AI have employed magnetic guides |15j | 
to limit transverse expansion of the sample and thereby 
extending the transit time. 
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FIG. 6. (Color online) (a) First order (N — 1) two-pulse echo signal for various applied B-field gradients, /3. The theoretically 
expected echo time is at At = t — 2Ti\ — 0. Each mA of current corresponds to a change of ~ 0.04 mG/cm in the applied 
gradient, (b) First order (N = 1) two-pulse echo energy as a function of /?. The solid red line is a fit based on Eq. (|12[) . The 
value of ft is obtained from a calibration using a flux-gate magnetometer. Pulse parameters for both (a) and (b): n = 800 ns, 
t 2 = 130 ns, I ~ 64 mW/cm 2 , T 2 i ~ 40.6 ms. 



B. Investigations of External £>-Gradients 

When T 2 i is large, the two-pulse AI can be used to 
explore the sensitivity to small external i3-gradients. We 
demonstrate the detection of changes in the B-gradient 
as small as ~ 4 x 10~ 5 G/cm in Fig. 6(a) 
N = 1 echo signal was recorded with T21 fixed at 



Here, the 
40.6 

ms for various applied gradients. Changes in the gradi- 
ent were facilitated by varying the current through the 
set of vertical quadrupole coils centered on the MOT (see 
Fig- El)- The smallest controllable increment in current we 
could achieve was 1 mA, which corresponds to a change 
of ~ 0.04 mG/cm as estimated from an independent cal- 
ibration based on a flux-gate magnetometer. 

In a similar experiment, the N = 1 echo energy was 
measured for T 2 \ fixed at ~ 40.6 ms as a function of j3, 
as shown in Fig. |6(b)| Here, it is clear that the echo 
energy has a strong periodic dependence on the applied 
B-gradient. These data provide confirmation of the the- 
oretical prediction given by Eqs. (0 and dHJ. This de- 
pendence is produced by the interference between elec- 
tric fields scattered off of gratings produced by different 
magnetic sub-levels. For example, for a given /3, gratings 
produced by states \Frrip) and \F m' F ) undergo phase 



shifts m,F<t>j^ 



and m'jr 



.(2) 
»/3 > 



respectively, where 



given by Eq. l[8|). For constructive interference between 

fields scattered by these states, the B-gradient must sat- 
(2) 

isfy (mF—Tn' F )(f>p /2 = 2n7r, where n is an integer. Thus, 
as (3 is varied, the phase shift induced in the \Fmp) and 
\Fm' F ) gratings produces periodic constructive (destruc- 
tive) interference in the total scattered field, and there- 
fore, maxima (minima) in the echo energy. This process 
occurs simultaneously in all 2F+1 sub- levels. As a result, 
the observed signal is a weighted sum of the scattered 
fields from all states. Here, there are 2F(2F + 1) = 42 



pairs of states that produce interference — although not 
all pairs have unique contributions. Since the excitation 
beams were circularly polarized in the experiment, the 
fields scattered from the extreme states (|3 3) or |3— 3)) 
dominate the signal. 

We use the following model, based on the squared mod- 
ulus of Eq. (J7J, to fit the data shown in Fig. |6(b)| 

S^(f3,T 21 )=S e-^-^ 2 / T2 a ™ F a m - F 

m F ,m' F (12) 
x e iA(m F -ra' F )l3N(N+l)(T 21 ~t 1 )' 2 

where So, ii and the set of {a mF } are free parameters, 
A = qgr/iB/ZM is a constant and to was set to T 2 \ for 
this data. In this model, the Gaussian factor outside the 
sum is added phenomenologically to account for signal 
loss due to both the transit time and any decoherence in 
the system. Also, each a mF is proportional to the mag- 
netic sub- level population, |a mF | 2 , through Eq. (|A21I) . 
As a result, these parameters are constrained to be pos- 
itive. All other fit parameters are unconstrained. In 
principle, it should be possible to obtain the sub-level 
populations from the set of best fit parameters {a mF }. 
However, determining the constant of proportionality be- 
tween the a mj7 , the populations and the scattered field 
intensity is complicated [H, H(| and not addressed by the 
theory presented here. We emphasize, however, that fits 
to data presented in this work give similar results for the 
set of {a mF }, which are consistent with our expectations 
for circularly polarized excitation beams. 

It is interesting that a measurement of the parameter 
A from data similar to that shown in Fig. 6(b)| can be 



used to test the theory of magnetic interactions [30J, [3 

Surveys of gradient-induced modulation on the echo 
signal shown in Fig. [7] provide additional confirmation of 
the theory outlined in Sec. [TT] and the Appendix. Fig- 
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FIG. 7. (Color online) (a) Two-pulse echo energy as a function of T21 in the presence of a fixed /3. Gradient-induced oscillations 
are shown for the first two echo orders (black points: N = 1; gray points: N — 2). Here, T21 is varied in integer multiples of 
T q to avoid sensitivity to atomic recoil effects. Fits based on Eq. (|12[) . shown as the solid lines, give |/3| = 9.30(1) mG/cm and 
/3| = 9.34(1) mG/cm, respectively, (b) Data analogous to (a) obtained using a slightly larger gradient with the three-pulse AI. 
Here, T21 is fixed at 2.0 ms and T32 is varied to map out the modulation for the first two echo orders (black points: N = 1; 
gray points: N = 2). Fits based on Eq. (fT3)l give \/3\ = 17.50(4) mG/cm and |/?| = 17.78(5) mG/cm, respectively. Two-pulse 
(three-pulse) AI parameters: t\ = 800 (800) ns; r 2 = 100 (70) ns; r 3 = (70) ns; / ~ 64 mW/cm 2 . 



ure 7(a) indicates that, in the presence of a -B-gradient, 
the two-pulse echo energy becomes modulated at a fre- 
quency that increases linearly with T21 (i.e. the modu- 
lation is chirped), as predicted by Eq. This figure 
shows gradient oscillations for both the N = 1 and the 
N = 2 orders of the two-pulse echo. Since the chirp rate 
increases as N(N + 1), the second order echo is modu- 
lated at a rate three times that of the first order echo. 
Confirmation of this is provided by a least-squares fit 
to the data based on Eq. (fl"2|) , as shown by the solid 
lines in Fig. 7(a) Since the gradient was held fixed in 



the experiment, the fits to the two data sets should pro- 
vide similar measurements of \(3\. The two measurements 
yield \/3\ = 9.30(1) mG/cm for N = 1 and \fi\ = 9.34(1) 
mG/cm for N — 2 [38], where the quoted error is the la 
statistical uncertainty generated by the fit. These mea- 
surements are in good agreement with each other and an 
independent measurement from a flux-gate magnetome- 
ter. We emphasize that accurate fits to these data and 
the extraction of j3 were possible only through the de- 
velopment of the multi-level formalism presented in the 
Appendix. In particular, since the oscillations shown in 



Fig. 7(a) do not occur with 100% contrast (i.e. each os- 



cillation minima does not reach the level of the noise) , a 
model including only two magnetic sub-levels with e qual 
excitation probabilities, such as that described in Ref . uft. 
is insufficient to model the data. 



Figure |7(b)| shows data similar to that shown in 
Fig. 7(a)| but for the first two orders of the three-pulse 
echo and a slightly larger .B-gradient. This data illus- 
trates that the three-pulse AI is less sensitive to gradi- 
ents than the two-pulse AI. Since T21 is fixed at 2.0 ms, 
the modulation frequency is constant and proportional to 
N and T21 — confirming the predictions of Eq. (ITT1) . The 



data is fit to the following model: 

^ 3 )(/?,T 32 ,T 21 ) = 5oe-( T -- t '') 2 / T2 Y, a ^ a m' F 

rriF ,7n' F 

xe iA( mF ~ m ' F )p[N(N+l)Ti 1 +2NT 21 (T 3 2-ti)]^ 



which is based on Eqs. (|A"25|) , (|A"26t and (|A"28]l . with 
a Gaussian decay factor added phenomenologically. All 
other parameters in this model are similar to those dis- 
cussed in reference to Eq. (|T2"j) . Measurements of the 
magnitude of the gradient from fits to these data yield 
\(3\ = 17.50(4) mG/cm and |/3| = 17.78(5) mG/cm for 
the N = 1 and N — 2 echoes, respectively. These two 
measurements differ by more than 5er, which deserves 
some explanation. By inspecting the fit to the N = 2 
echo, it is clear that the data is not well-modeled by a 
single frequency sinusoid as T 32 becomes large. This pro- 
vides evidence that the atoms are sampling different gra- 
dients as they drop under gravity — an effect that is not 
accounted for in the theory. By analyzing different sec- 
tions of this data, we estimate that the gradient varies by 
as much as ~ 1.6 mG/cm between T32 ~ 40 ms and 100 
ms — during which time atoms fall ~ 4 cm. Independent 
measurements of the curvature of the -B-field, where \j3\ 
was found to change by ~ 0.4 mG/cm every centimeter, 
are consistent with the variation in j3 detected by atoms. 

Although we have demonstrated sensitivity to changes 
in the -B-gradient as small as ~ 4 x 10~ 5 G/cm using 
T21 ^40 ms with the two-pulse AI, our ability to mea- 
sure the absolute magnitude of the applied gradient is 
less sensitive. This is primarily because the measure- 
ment is based on fitting data to an oscillatory model and 
extracting the modulation rate — which cannot be done 
accurately without the presence of an oscillatory com- 
ponent in the data. To estimate the smallest measur- 
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FIG. 8. (Color online) Smallest direct measurement of a 
B-gradient using the two-pulse (a) and three-pulse (b) AIs. 
Here, the applied fields were tuned separately for each AI 
such that the first revival in the N = 1 echo energy occurred 
at the largest time. A fit based on Eq. (|12|l for the two- 
pulse AI gave |/3| = 0.26(3) mG/cm. Similarly, a fit based on 
Eq. (JT3J yielded |/3| = 9.5(1) mG/cm for the three-pulse AI, 
with T21 ~ 1.27 ms. Two-pulse (three-pulse) AI parameters: 
Ti = 800 (800) ns, r 2 = 100 (70) ns, r 3 = (70) ns, / ~ 64 
mW/cm 2 . 



able B-gradient with the two interferometers, we tuned 
the applied fields for each AI separately such that the 
first revival in the N = 1 echo energy occurred at the 
largest time. The resulting data are shown in Fig. |8j 
which yielded measurements of \f3\ = 0.26(3) mG/cm for 
the two-pulse AI [Fig. 8(a)] a nd 
the three-pulse AI [Fig. |8(b)| . 



= 9.5(1) mG/cm for 



V. APPLICATIONS TO GRAVITY 

The apparatus shown in Fig. |3j is designed for mea- 
surements of the atomic recoil frequency 26]. As a re- 
sult, it is not isolated from external vibrations and is 
unsuitable for measurements of the optical phase of the 
scattered read-out light using heterodyne detection. For 
this reason, a measurement of g from the phase of the 
atomic grating [22], (27j is beyond the scope of this Article 
and will be presented elsewhere. However, the aforemen- 



tioned results relating to £?-gradients validate theoretical 
predictions that can be applied to precise measurements 
of gravity. In this section, we discuss the feasibility of 
such a measurement by applying the formalism presented 
in the Appendix. 

The best portable gravimeter [3{| uses an optical 
Mach-Zehnder interferometer where one arm contains 
a free-falling corner-cube for position-sensitive measure- 
ments of g at the level of ~ 1 ppb over a few minutes. The 
position sensitivity in these devices comes from detecting 
interference fringes as a function of the drop time of the 
cube relative to an inertial frame defined by a stationary 
mirror. The frequency at which the fringes accumulate 
scales linearly with the drop time (i.e. the frequency 
is chirped). The matter- wave analog of this gravimeter 
is the two-pulse echo AI [22], [27j , where changes in the 
phase of the grating due to gravity are detected relative 
to the nodes of a pulsed sw — which serves as the inertial 
reference frame. In this case, the accumulation of fringes 
due to matter-wave interference is also described by a 
chirped-frequency sinusoid. 

We now review the main results of the grating echo 
theory that pertain to gravity. The gravitational poten- 
tial can be written as 



U{z) = Mglz, 



(14) 



where the force is T = — Mg and I is the (2F + 1) x 
(2F + 1) identity matrix. The effect on the echo AI is 
similar to that of the _B-gradient on a sample that has 
been optically pumped into a single state. Since gravity 
acts equally on all states, the phase shift of the grating 
produced by each state is the same. Therefore, the ex- 
pression for the field scattered from the grating simplifies 
significantly compared to Eq. 



Ef\t;T) 



£4?l(*;T) 



(15) 



where the grating phase due to gravity is 



\At;T 21 ) = -f{N(N+l)Tl 
+ 2 [Tr + (N + 1)T 21 ] At + Ai 2 } 



(16) 



as determined by Eq. (|A23I) . This phase cannot be de- 
tected from the intensity of the scattered light because 
there is no differential phase shift between magnetic sub- 
levels — thus, there is no amplitude modulation of the 
grating [27| ■ The scaling of the grating phase with T 21 in 



Eq. ([TB"]) shows the similarity between the two-pulse AI 
and the optical Mach-Zehnder interferometer discussed 
in Ref . M 

Figures 9(a) and |9(b)] show the expected two-pulse AI 
signal in the presence of gravity as a function of T%\ — 
illustrating that the modulation frequency is chirped lin- 
early with T 2 i (<4 2) = d(j} { g ] /dT 21 oc T 21 ). As T 2 i in- 

(2) 

creases, uj q becomes larger than the recoil frequency (for 
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FIG. 9. (Color online) The real part of the scattered electric 
field (red solid line) in the presence of gravity for the two-pulse 
(a,b) and the three-pulse (c-f) AIs. In parts (a-d), the black 
dashed line shows the signal envelope [based on Eq. (|A2I|) 
for (a,b) and Eq. (|A29[) for (c,d)], which is modulated at 
2ui q ~ 2n x 31 kHz. (a) The field oscillates at a frequency 
w ( g 2) (T 2 i) ~ 2tt x 50 kHz in the vicinity of T 2 i ~ 1 ms. (b) 

(2) 

The modulation frequency increases to tjg ~ 27T X 450 kHz 
at T21 ~ 9 ms. Parts (c,d) show results for the three-pulse 
AI as a function of T21, with T32 fixed at 100 /xs. The field 
oscillates at approximately the same frequency in (a,c), and 
in (b,d), since T32 is small. Parts (e,f) show the three-pulse 
signal as a function of T32 with T21 = 5 ms. Since T21 is 
fixed, there is no sensitivity to atomic recoil and the signal 
envelope is not modulated. As T32 increases, the modulation 
frequency of the scattered field amplitude remains fixed at 



,(3) 



(T2 



2tt x 125 kHz. 



tion frequency is given by 

= qg[N(N + 1)T 21 + NT 32 + (N + l)At], (17) 



4 3) 



dT 21 



which is identical to the two-pulse case, 



.(2) 



with the 



addition of the term proportional to T32. 

Figures 9(e) and |9(f)| show the expected three-pulse 
signal as a function of T32, with T 2 \ fixed at 5 ms. In 
this case, there is no sensitivity to atomic recoil, so the 
envelope remains at a constant level as T32 is varied. The 
frequency of the phase modulation is also fixed by N and 
T 2 i, as given by 



ui 



(3) 



(T 2 



8T- 



:S2 



qg (NT 21 + At) 



,(3) 



(18) 



2ttx 



For the conditions presented in these figures, U) g 
125 kHz. The work of Ref. UH shows that the three- 
pulse AI is significantly less sensitive to mirror vibrations 
than the two-pulse AI if T32 ^ T 2 \. Our results have 
also shown that this configuration is less sensitive to B- 
gradients. For all these reasons, this AI is particularly 
well-suited for precise measurements of g. 

Simulations of the two-pulse AI signal with N = 1, 
T21 ~ 150 ms and a phase error of 1% suggest the preci- 
sion of a measurement of g should be ~ 1.4 ppb. Simi- 
larly, we estimate a precision of ~ 0.4 ppb for the three- 
pulse AI using N = 1, T 2 i = 75 ms, T 32 varied over 150 
ms and the same phase error. From these estimates, it is 
clear that these AIs can have greater sensitivity than the 
best industrial sensor [39(. Since the precision scales lin- 
early with the phase error, we anticipate further improve- 
ments in sensitivity without extending the timescale. If 
systematic effects of such a cold atom gravimeter are 
characterized, it may be possible for the AI experiment 
to serve as a reference to calibrate other gravimeters. 



the first order echo in 85 Rb, this occurs when T 2 \ > 300 
/is), and T21 must be incremented in steps less than r q 
to avoid undersampling the frequency. However, this 
effect causes reduced sensitivity to the grating phase, 
since modulation at the recoil frequency produces pe- 
riodic regions with small signal amplitude. Additionally, 
as shown in Ref. [HI, this AI is very sensitive to phase 
changes due to mirror vibrations, which can be detri- 
mental to measurements of g using this technique. 

Figures 9(c) and |9(d)| show the expected three-pulse 
signal as a function of T 2 \ in the presence of gravity. 
It is clear that the envelope of the scattered field has 
a complicated periodic dependence on T21, with a zero 
every r q ~ 32 /js due to the destructive interference of 
momentum states differing by the two-photon recoil mo- 
mentum, hq. This is similar to the two-pulse case shown 
in Figs. 9(a) and |9(b)| Here, the grating phase modula- 



VI. CONCLUSIONS 

Measurements of applied i3-gradients using both the 
two- and three-pulse techniques arc in good agreement 
with independent measurements of ft using a flux gate 
magnetometer. We have demonstrated sensitivity to 
changes in the i?-gradient at the level of ~ 4x 10~ 5 G/cm. 
Absolute measurements of ft as small as ~ 3 x 10~ 4 G/cm 
were also possible using the two-pulse AI. These measure- 
ments indicate that an accurate description of the data 
presented above requires the inclusion of multiple mag- 
netic sub-levels. We have also shown sensitivity to spatial 
variation in the B-gradient using a long-lived second or- 
der (N = 2) three-pulse echo. It is this non-linearity in 
the -B-field that affects the timescale in echo AIs rather 
than the presence of a small, uniform B-gradient. 

As tests of the theoretical results presented in Sec. [Til 
we have separately confirmed the linear dependence of 
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the /3-induced oscillation frequencies, ' and uj^ ' [given 
by Eqs. ((9|) and (fTTj) . respectively], on the £?-gradient. 
We have also verified that these frequencies both scale 

(3) 

linearly with T21, and, for the three-pulse AI, is con- 
stant as a function of T32. 

Since we have achieved signal lifetimes approaching the 
transit time limit, we have shown that fountain-based 
experiments are possible with grating echo AIs. The ad- 
vantage of a fountain configuration is that the spatial 
extent of the AI (~ 11 cm for 300 ms timescale) can 
be made small, which reduces the requirements for in- 
homogeneous B-field suppression. Such a configuration 
is ideal for precise measurements of gravity, particularly 
with the three-pulse AI. Passive suppression of -B-fields 
with larger cancelation coils, or optically pumping into 
the rriF = sub-level, represent two ways in which such 
a measurement can be realized. 

Despite the widespread use of Raman-type AIs for in- 
ertial sensing H,|13,|4l|, grating echo- type AIs — which of- 
fer reduced experimental complexity — are also excellent 
candidates for precision measurements of to q and g. This 
work has brought about understanding of systematic ef- 
fects produced by -B-gradients on these measurements. 

In summary, we have developed a complete under- 
standing of the effects of a constant force that applies 
to all time-domain AIs. Although the sensitivity for AI- 
based gradient detection cannot compete with commer- 
cial magnetic gradiometers (which offer sensitivities of 
the order of ~ 1 pT/m), the technique is useful for abso- 
lute measurements of gradients in cold atom experiments. 
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APPENDIX 

In this appendix, we derive expressions for the signals 
generated by the two- and three-pulse interferometers in 
the presence of a constant external force, J- . In Ref. [27l . 
a similar calculation for the two-pulse signal is given, 
in which only two ground state sub-levels are consid- 
ered, and effects due to spontaneous emission are ignored. 
Here, we account for 2F + 1 magnetic sub-levels in the 
field scattered from the atoms, as well as spontaneous 
emission during the excitation pulses. Both of these ef- 
fects are crucial for an accurate description of these in- 
terferometers. We also give a general expression for the 
signal generated by an TV-pulse AI from which all classes 
of grating-echo interferometers can be realized. 



The potential is assumed to have the form U(z) = 
—Mz, where M = —dU/dz is an operator that computes 
with z and p, and acts on the basis states \F vtif) with 
eigenvalues mpj- '. Here, J 7 is a constant with units of 
force. We proceed by computing the ground state wave 
function after the application of each sw pulse at times 
t = T\ and T2, with a period of evolution before, between 
and after each pulse (with durations T\, T2 — T\ and 
t — Tii respectively) in the presence of the force. During 
the application of each sw pulse, the kinetic and potential 
energy terms in the Hamiltonian are ignored by assuming 
the pulses are sufficiently short such that the atom does 
not move significantly (Raman-Nath approximation). In 
this manner, the sw pulses are treated as Dirac (S-function 
excitations, even though they are given durations Tj for 
the purposes of the calculation. 

The interferometer signal is defined as the back- 
scattered electric field amplitude at the time of an echo, 
which is proportional to the amplitude of the q-Fouricr 
harmonic of the density distribution at these times. The 
results for the two-pulse AI signal are then generalized for 
an A^-pulse AI, from which we compute the three-pulse 
AI signal. 

The Hamiltonian for the ground state \F mj?) in the 
presence of a sw field and an external potential, U[z), 
can be approximated by [TH, [ItJ 

2 

H mF = 7^ + f>-Xm F e ie cos(qz) + U(z), (AI) 

where 9 is a phase associated with spontaneous emission 
during the sw pulse 

9 = tan" 1 (-■!), (A2) 

and Xm F is a two-photon Rabi frequency given by 

X™F = 2A \ + J (PmP 1L m F +q L ) ■ (A3) 

Here, fl is the on-resonance Rabi frequency for a two- 
level atom, A = u>l — loq is the atom- field detuning 
with atomic resonance frequency ujq and laser frequency 
u>l, 7 is half of the spontaneous emission rate, and 
(Cm F q L mt+qt) ^ s a Clebsch-Gordan coefficient for a 
light field with a polarization state qL- We ignore the 
excited state in this treatment, since the field is assumed 
to be relatively weak and far off- resonance (| A| 3> Qo, 7). 
We also neglect the Zeeman shift of magnetic sub-levels 
by assuming |A| 3> gp^sB/h. 

The amplitude of the ground state wave function at 
t = can be written as a superposition of spin states: 

0) \F m F ) , (A4) 

mp 

where the amplitude of each spin state is 

a mF (z,0) = ^=e lp ° z /\ (A5a) 
a mF {p, 0) = a mp 5(p-p ). (A5b) 
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Here, po is the initial momentum of the atom along the z- 
direction, \a mp | 2 is the population of state \F nip), with 
J2m F \ a m F \ 2 — 1, and a mF {p, 0) is the amplitude of the 
spin state in momentum space. 

The main challenge in this calculation is evolving the 
wave function between sw pulses in the presence of 
the additional potential energy, U(z). In the absence 
of this potential, it is straightforward to integrate the 
Schrodinger equation in momentum space. However, 
with U (z) present, we have the following equation of mo- 
tion: 



e A e B , but one can use the Zassenhaus formula [HJ to 
expand the matrix exponential as 



MA+B) 



e ZA e ZB e -^[A,B]/2 



x e 



} ([A,[A,B]]-2[[A,B],B])/6 .. 



(A8) 



Mz a mF (p,t). 



dt V 2M 
One can integrate this equation to find 



(A6) 



(A7) 



where £ is an arbitrary constant. The higher order fac- 
tors (represented by • • • in the above equation) vanish if 
[[A, B],B] and [A, [A, &]] commute with all higher or- 
der nested commutators. Choosing A = —Mz and 
B = p 2 /2M [Hj], and using the commutation relations 
[z,p 2 ] = i2hp, [z,p] — ih, we find: 



-Mz, 



Mz, 
2M'' 



P 

2M 

- V 2 
-Mz,^— 
' 2M 



-Mz 



P 

2M 



but some care must be taken when evaluating the right 
hand side. The challenge arises from the fact that z and 
p = —ihd/dz are non-commuting operators. As a result, 
the exponential in Eq. (|A7[) is really a matrix exponential 

of non-commuting matrices A and B. In general e A+B 7^ 

I 

Using Eq. (|A8[) with £ = —it/K and the commutators in Eqs. (|A9I) . Eq. (| A7|) becomes 



hM 

h 2 M 2 



= 0. 



(A9a) 
(A9b) 
(A9c) 



(A10) 



Since e^'^' \F mp) — e^" lFjr ^ \Fmp), it follows that the amplitude of the state \F mp) before the onset of the 
first sw pulse is 



( Pi t) = amFe *(™FJ r )tz/h e -ip 2 t/2Mh e -i( mF r)pt 2 /2Mh e -i( mF r) 2 t 3 /6Mh S ^ p _ po ) 5 

rriF 

2tt% 



a m (z t) = QmF p i(Po+mFJ r t)z/h p -i<iot/h p -Hm F r)pot 2 /2Mh p -i(m F F) 2 t 3 /6MH 



(Alia) 
(Allb) 



where e = p\j2M is the initial kinetic energy of the atom. 

The first sw pulse, applied at t = T\, diffracts the atom into a superposition of momentum states. The wave function 
is computed in position space using the Raman-Nath approximation and integrating the Schrodinger equation to obtain 



ag\{z,T x ) = a mF (z,T 1 )^(- J )"J„(9W )e m « 



tt^ ) (p,T X ) = amFe -^oT 1 /h e -t(,yi F ^)p T 2 /2Mh e -i( mF r) 2 T 3 /6Mh 



(A12a) 
(A12b) 



(-i) n J n (0W )<J(p - po - mpFTi - nhq). 



Here, Qm F = «Eie' e is the (complex) area of pulse 1, u 



(i) _ 



,,,, — Xm F T i, T i is the duration of the pulse, and a$ F (p, T\) 
is the wave function in momentum space. The superscript (1) on a,m F denotes the number of sw pulses that have 
been applied to the atom so far. We use the prescription of Eq. (|A10[) to evolve the amplitude in momentum space 
[Eq. (|A12b|) ] until the onset of the second pulse 



(p )t ) = amFe »(m F ^)(t-T 1 )z/R e -i[p2T 1 +p 2 (t-T 1 )]/2M?i e -i(m F ^)[poT 1 2 +p(t-T 1 ) 2 ]/2Afr l 



(A13) 
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To apply the next sw pulse to the wave function, it is convenient to transform back to position space: 

<2 (1) (z t) = a "" lF c i (PQ+ m FFt)z/h c -i\p?,T 1 + (p ~ mF rT 1 ) 2 (t~T 1 )]/2Mh 



. 2nh 

x e -i( mF r)[p T? + (p +rn F FT 1 )(t-T 1 ) 2 ]/2Mh e -i(rn F r) 2 [T? + (t-T 1 ) 3 ]/6Mh (A14) 
x ^(_j)«J n (ea) ) e m9z e -mg«o(t-T 1 ) e -m 2 tl ;,(t-T 1 ) e - 4 n ? (m F ^)[(t-T 1 ) 2 +2T 1 (t-T 1 )]/2Af^ 
n 

Here, vq — Po/M is the initial velocity of the atom and ui q — hq 2 /2M is the two-photon recoil frequency. Applying 
the second pulse at t = T2, the wave function becomes 

n( 2 ) (y T \ — r i(Po+rn F TT2)z/h-i£ 2 ,T2/h r -ipo(m F r)T 2 /2Mh r -i(m F F) 2 Ti/6Mh 
a m F V^j ± 2) — r- — - e e e e 

y/ZTTn 

x yv_f)(«+™) j„(ew )j m (e, (2) ) e i ("+ m )9 z e ' m? "» (r! - Ti V i " 2u « (T2 " Ti V i "'' (mF - F)(T22 " T i )/2M . ( ' A15 ' > 

To evolve the wave function in the presence of the external force until time i, once again we transform into p-space 
and use Eq. (|A10p to obtain 



a $ (p,t) = amFe Krn F F)(t-T 2 )z/h e ~i[p 2 T 2 + P 2 (t~T 2 )]/2Mh e -^ 

x ^ H) (n+ m ) Jn(e M )Jm(e Wj e -^ (A16 ) 

x S[p — Po — mpTT-z — (n + m)Hq]. 
Finally, the amplitude in position space after the second pulse can be shown to be 

(z t) = a " lF c i (Po+ m F^t)z/h c -ieot/h c -i(m F J^)pot 2 /2Mh c ~i(m F T) 2 t 3 /6Mh 

mF \f2~nh 

X^(-z)("+ m )j4el 1 ^)J m (e^) F ) e ^«+™)^ e ~^o[™(T 2 -T 1 ) + (n +m )(t-T 2) ] (A1?) 



n . in 



x e 



-iw ? [n 2 (T 2 -Ti) + (n+m) 2 (i-T 2 )] -iq(m F F)[n(T 2 -T 2 ) + (n+m){t 2 -T 2 )]/2M 



To compute the field scattered from the atomic interference as a function of t, we use the q- Fourier component of the 
ground state density, Pm F m F {z,t) — \a,m F (z,t)\ 2 , which can be shown to be 

71,771,71/ ,771' 

x e -i< ? «o[(ri-nO(r 2 -T 1 ) + («+m-n'-m')(t-T 2 )] e -it J ,{(r l 2 -n' 2 )(T 2 -T 1 ) + [(n+m) 2 -(n'+ m ') 2 ](t-T 2 )} ( Al8 ) 
x e -iq(m F r)l(n~7i')(T 2 2 ~T?) + (n+m-n'-m')(t 2 -T 2 2 )]/2M 

Since the density distribution contains frequency components that depend only on the difference between interfering 
momentum states, we recast the sums over n' and m' in terms of vN = n — n' and v = nl + m 1 — n — m (the integer 
difference between momentum states after the first and second pulses, respectively): 

PijiFraF^i'i " 2 n fi 2—< J nV<J mF ) J n-vN\ KJ rn F ) J m\ VJ m F ) J m+v(N+l)\ KJ m F ) e 

x e iuqv [(t-T 2 )-N(T 2 -T 1 )] e iuu q {[2(n+m)+u](t-T 2 )-N(2n-uN)(T 2 -T 1 )} ^ > 

x e iuq(m F F)[(t 2 -T 2 )-N(T 2 -T 2 )]/2M _ 

The scattered field is proportional to the g-Fourier harmonic of jOmlfmr^, i) [the coefficient of the e~ luqz term in 
Eq. ([Mil) , with y = !]• 

Summing over all magnetic sub-levels in the ground state, one can show that 

E% ] (t; T) = E ^ml (*; T)e lmF ^ ^ , (A20) 
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where 



iV 

x (2«« ^sin^ - 0) sin(¥>i + »)) (2«<& Vsin(^ 2 - 0) sin(^ 2 + 0)) (A21) 

'8in(^+e)\^Vsin(^-e)\^ +1)/2 



v sin(<^i — 6>) / \sin(y> 2 + 0) , 
is the field scattered from each magnetic sub-level, with recoil phases 

Mt;T)=u; q {t-t^l ), (A22a) 
<f 2 {t;T) =u q (t-T 2 ), (A22b) 

(2) 

and m,F<pj: is the phase shift of the density grating produced in the ground state \Fitif) due to the presence of the 

(2) 



external force, J 7 , with cjy^' given by 

^ (<; T) = ui [(<2 " Tl) " * (T » " Tl2)] ■ (A23) 

In deriving Eq. (IA21[) we have made use of the Bessel function summation theorem 

HE Eft m 

£ Jn {ue lf) ) J n+V (ue- ie ) e « 2n +^ = iV„ (2^^/^ - 0) sin(0 + 0)) ^^1 ^ , (A24) 

and we averaged over the velocity distribution of the sample assuming a Maxwellian distribution centered at vq with 
e _1 width <j v = y / 2ksT /M . In this way, we account for the possibility of an initial launch of the atomic cloud 
and for the dephasing of the echo due to the distribution of Doppler phases in the sample. An additional factor 

mp+,J w as added to the scattered field to account for the atom-field coupling by the read-out pulse. 
The scattered field lasts for a time r co h = 2/qa v — called the coherence time — about each echo, which occur at times 
^icho = N(T2 — Ti) + T2. The phase in Eq. (|A21I) . associated with spontaneous emission during the excitation pulses, 
affects only the recoil-dependent component of the signal [ljj ■ 

These results can be generalized for the case of an A-pulse interferometer with a set of onset times T = 
{Ii, T 2 , ■ ■ ■ , T N } for which T J+1 > Tj. After N sw pulses, each with pulse the total scattered field at 

time t is 

E ( p (t; T) = £ E^J (*; T)e im '+™ (t;T) (A25) 

where 

' M<Pi-e) \ lh - l '- l)/2 iA2(i! 

sm(cpj + 6) 



II J (h-h-i) ( 2u rn F \/ sin (^ - 8) sinifj + 
3=1 



Here, I = {l±, 1%, . . . , In} denotes the set of momentum states that interfere after the pulse sequence, where lj is the 
difference between interfering momentum states (in units of hq) after pulse j. The echo times and the recoil phases 
are given by 



N-l 

N 

<pj(t; T) = uj q h{T k+ i - T fc ), (A27b) 
k=j 
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and the contribution to the phase of the grating due to the force, T, is 

N 

^(^)^EM^^ 2 )- (A28) 
j=i 

In Eqs. (|A26[) - (|A28I) Ijy = 1, which corresponds to the scattered field from the q- Fourier harmonic of the density 
formed after the sw pulses, and it is understood that lo = and T^+i = t. 

We now use the formalism for the iV-pulse echo signal [Eq. (|A26[) ] to obtain an expression for the three-pulse 
interferometer signal discussed in Sec. HU We begin by setting N = 3 and T = {Ti,Ti + T 2 i,T 1 + T 21 + T 32 }. For 

an echo_to occur at t { c % = T x + T 32 + (N + l)T 2 i for any Ty, T 32 and T 21 , Eq. (|A27ap dictates the set of lj to be 
I = {— N, 0,1}. Then, it can be shown that the scattered field is given by 

N 

x J R (2u£l Vsin^i - 0) sin^i + 0)) J R (2u<& ^sin& 2 - 9) riKfo + 0)) (A29) 



x Ji (2u££ Vsin(^3 - 0) sin(^ 3 + (9) 
where the recoil phases in this case are 



sim>! + 6) \ R ' 2 ( sin(v3 2 -6)\ R/2 ( sin^a - ^ ^ 1/2 



sin^i — 0) / \sin(( y 9 2 + Vsi 11 ^ + 



and the grating phase due to T is 



<Pi = w, (t - ifl) , (A30a) 
^2 = <P3 = ^ g (< - t£> + 7VT 21 ) , (A30b) 



= ^ {N(N + l)Tl + 2NT 32 T 21 + 2[Ti + T 32 + (N + 1)T 21 ] At + At 2 } . (A31) 
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